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We present an internal shocks model to investigate particle acceleration and radiation production 
in microquasar jets. The jet is modelled with discrete ejecta at various time intervals. These ejecta 
(or 'shells') may have different properties including the bulk velocity. Faster shells can catch up 
and collide with the slower ones, thus giving rise to shocks. The particles are accelerated inside 
the shocked plasma. Each collision results in a new shell, which may take part in any subsequent 
collisions as well as radiate due to synchrotron radiation. Almost continuous energy dissipation 
along the jet can be obtained with a large number of shell collisions. We investigate the spectral 
energy distribution of such jets as well as the physical significance of various parameters (e.g. the 
time interval between ejections and the shell size) 
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1. Introduction 



Jets are observed from a large variety of astronomical objects, ranging from young stellar 
objects (YSO) to active galactic nuclei (AGN). In terms of emission, flat radio (a = when Fy o<: 
v") spectra have commonly been observed from these sources (AGN and X-ray binaries (XRB)). 
Observations of XRBs suggest a frequency independent flat spectrum extending down to millimetre 
[Q] or even near-infrared bands |j3|]. Flat radio spectra from AGN and XRBs have traditionally been 
attributed to partially self-absorbed jet emission from the innermost region of a conical jet [jl|], 
However, in order to obtain such flat spectra from partially self absorbed synchrotron emission 
we either need continuous re-acceleration all along the emitting jet region or use a non conical jet 
model If we assume a conical geometry for jets then the acceleration mechanisms problem 
needs to be addressed. We explore internal shocks as a re-acceleration mechanism for microquasar 
jets. 



2. The model 

Our model has been inspired by the internal shocks model for radio-loud quasars [^. We 
assume a conical geometry for the jet with a given full opening angle dg- Plasma is injected into 
the jet in the form of discrete packets or "shells". 

2.1 Shell properties 

Our model is heavily based on Spada et al. model ^ with a few changes in order to adapt the 
model to XRBs. The injected shell properties determine most of the observables for the simulation. 
The central engine for shell ejections is considered to be intermittent in nature, with a time interval 
(A,) between ejections. Each shell also has a given bulk Lorentz factor (BLF), F, associated with 
it. Outside of the collisions, the F does not evolve as a shell moves through the jet. In order to 
choose F for a given shell we set the minimum(F„„>, > 1) and maximum(Fma.v) T and then sample 
randomly from that range. 

The shell mass is determined by the combination of kinetic luminosity of the jet {L\f), the time 
jet is on for {tjet) and the total number of shells to be injected during that time (A'^). The following 
relation sets the total mass in the jet: 

N 

Y,MjrjC^=Lwtjet , (2.1) 

7=1 

where c is the velocity of Ught and Mj is the mass of an arbitrary shell. However, as we have 
assumed the central ejection engine to be intermittent in nature and accumulative, the time gap 
between ejections simply accumulates more mass for a given shell. 



Mj = ^, (2.2) 



L^r/\t j 

where Atj denotes the time interval between ejections j — I and j. We also set the width of the 
shefl (Ay). 
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In order to calculate the time gap between ejections, we take a Gaussian distribution with the 
mean corresponding to the mean time gap ( Afy) we desire. We then sample from that distribution 
with chosen standard deviation. If the standard deviation is set to be low then the average time gap 
between ejections will be close to the mean and the total number of shells to be injected will be 
« fjef/Af. With a larger standard deviation, the total number shells to injected will deviate from 
the previously stated relation. The variations in the ejection time gap as well as the bulk Lorentz 
factors of the shells leads to various shells catching up with each other, causing collision. 

2.2 Two-shell collisions 

Once the above parameters are determined, we can set the injection times as well as all the 
properties of the shells at the start of the simulation. With the variations in T and /\tj we can have 
a large number of collisions taking place. The time gap between any two collision is: 

ninner pouter 

^ [A;-i)-i8# + 0-5[i8(%i)+i8/]c ' ^^"^^ 

where R'""'^'' and R"""''' are the radii to shell inner (with respect to the source) and outer boundary 
i.e. the shell has a finite width. j8 is the shell velocity and jS*^ is the shell thermal expansion velocity, 

R<' = ^ I (2 4) 

where jS^c = v'^ is velocity of sound in the plasma in the shell co-moving frame with 



Eth denotes the thermal energy of the shell, causing the shell to expand in between collisions. 

The collisions are treated as inelastic collision. Momentum and energy conservation laws give 
us the merged shell properties. Merged shell mass is simply given by: 

M,„=M,+Mo, (2.6) 

where subscripts / and o denote inner faster and outer slower shell involved in the collision. If we 
have the interacting shells bulk Lorentz factor. Mass and internal energy (t]), the merged shell bulk 
Lorentz factor has the form: 



^At,-/r,-+At«/r« 

where jXi = Mi + rji/c^ and /Iq = Mo + f]o/c^- The internal energy of the merged shell is given by: 

Ein = rii + Tlo + At/c^(r,- - r^) + ^oC^iTo - r„). (2.8) 

The merged shell width is calculated by hydrodynamical treatment of the shocked plasma. Adapt- 
ing the jump equations [Q] for forward and reverse shock at the point of collision of the two shells 
we get the width (in the lab frame) of the outer shock as: 
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A.,,,o = ^ , (2.9) 

Po 

where Ao is the outer shell width and 

_r,„ yr, + i 

Po-p^ 7-1 ' ^^-^"^ 

where 7 is the adiabatic index. Equivalent equations stand for the reverse shock involving the inner 
shell. The width of the two corresponding shocks then gives the width of the newly formed shell. It 
should be noted that the overall effect of the above treatment is to cause compression in the plasma 
involved in the collision. 



2.3 Shocked plasma 

The energy generated in the shocked compressed plasma is used in the following three ways 
(assuming equipartition): 



1. Accelerate electrons to a powerlaw distribution. 

2. Generate magnetic field for the synchrotron radiation. 



3. Thermal expansion of the shells. However, we do not explicitly take thermal electron distri- 
bution into account. 



The energy split between the above is parameterized in the form of Ue, ub and Uth- These determine 
the fraction of the total energy generated that can be made available for the respective process (i.e 
Ue + UB + Uth = Eint)- The parameter Uth determines the thermal energy content of a given shell 
which in turn determines the local sound speed. This is used to calculate the longitudinal (along 
the jet) expansion velocity of a shell. 



2.4 Powerlaw electron distribution 

With each collision, the electrons are accelerated to a powerlaw distribution of the form: 



N{E)dE = KE-PdE , (2.11) 

where E = yrUeC^ is the electron energy and K is the normalization constant (which can be obtained 
by integrating the above relation between Emin and Emm:)- The acceleration is treated as instan- 
taneous and uniform in the respective volume. At each collision even the powerlaw distribution 
of electrons is completely replaced by a new one. However, in between collisions the powerlaw 
distribution "ages" via the evolution of K as well as Emax- 



With the assumption of a conical jet and expanding shells (see 2.2), we need to take adiabatic 



losses into account as well (see 2.6). 
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2.5 Magnetic field 

Parameter ub determines the fraction of energy available for the magnetic energy (Eb = us-Eint): 

Eb = ^: (2.12) 
2iUo 

where /Iq is the magnetic permeability. This gives the magnetic field strength available to syn- 
chrotron radiation. The magnetic field is assumed to be randomly oriented and tangled in the 
plasma. This allows us to treat the magnetic field as an ultra-relativistic gas (7 = 4/3), giving a 
straightforward adiabatic loss treatment for the magnetic field. 

2.6 Adiabatic losses 

As we are dealing with expanding shells, we need to take the adiabatic losses into account. The 
shells expand both laterally due to conical jet as well as longitudinally due to thermal expansion 
of the shells. Using the treatment outhned in [^, we get a relation for the change in energy due to 
adiabatic losses 

E fV\-'''^ 

(2.13) 



£0 \Vo 

This relation can therefore be used for both electron energy losses as well as the magnetic energy 
losses. 

2.7 Synchrotron spectrum 

The Synchrotron radiation treatment outlined in is employed in our model. As noted earlier, 
we assume that any given shell moves with constant velocity along the jet, unless it collides. For a 
shell with bulk Lorentz factor T and velocity j3 and the jet axis making an angle d with the observer, 
we get the Doppler factor of 

5^ = [T{\^^cose)]-' , (2.14) 

where =F are the receding and approaching jets respectively. This then leads to observable monochro- 
matic intensity 

where Jy is the emissivity per unit volume, Xv is the absorption coefficient and r is the radius of the 
shell. The frequency is measured in the rest frame of the jet material. It is related to the observed 



frequency by Vohs = 5^v. With the estimate of distance to the source and equation 2.15 , we can 
calculate the flux density for any given shell. 

3. Results 

As mentioned earlier, the Radio to IR flat spectrum from a conical jet is a commonly accepted 
picture (for XRBs and AGN), yet not fully explained. Therefore, with our model in its early 
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Stages of development, obtaining the flat spectrum is quite an interesting diagnostic. However, it 
is highly non-trivial to obtain a realistic flat spectrum from a conical jet without making certain 
assumptions about the re-acceleration mechanism. We present some preliminary results that show 
the self consistently calculated synchrotron spectra obtained from a conical jet with internal shocks 
taking place along the entire length of the jet. 
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Figure 1: Radio and infra-red lightcurve from the entire jet. The average fluxes and spectral index are shown 
with straight lines. 



Figure || shows the radio and infra-red lightcurve from the entire jet. The lightcurves not only 
shows how the infra-red and radio flux vary with time, but also the spectral index. It can be seen 
that the spectral index oscillates between a > (inverted) and a < 0. The spectral index can also 
be seen to spend time in the flat spectrum (a = 0) regime. Figure ^ shows the spectrum from jet 
a with almost identical parameters as used for figure [T], higher kinetic luminosity (the parameters 
are summarized in table [T]). It can be seen in figure ^ that IR flux is much higher than radio, giving 
rise to an overall inverted spectrum. The possible reasons are discussed in the next section. For 
completeness figure |3| shows a lightcurve from a weak jet. This is primarily to show that with low 
energy densities, the shell are almost instantly optically thin to radio and hence the radio flux is 
greater than the IR flux. 

The key parameters giving rise to the spectra are summarized in table |l[ The injection fre- 
quency is the average time gap between shell injections and y,nax is the integration upper limit for 



powerlaw electron distribution (equation 2.1 1 ). Similarly T„^x is the upper limit for BLF sampling 
range. 
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Figure 2: Spectrum showing IR flux much greater than the Radio flux. Average flux and spectral index 
values are also shown. 
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Figure 3: A lightcurve form a weak jet. Average flux and spectral index values are also shown. 
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Table 1: Model parameters to achieve the spectra shown (see Results section for details) 



Parameters 


Figure 1 


Figure 2 


Figure 3 


Jet kinetic luminosity (erg s^^) 


1 X 10^^ 


1 X 10^' 


1 X 10^^ 


Distance (kpc) 


2 


2 


2 


Injection freq. (Hz) 


1 


1 


1 


Ub 


0.1 


0.1 


0.1 


Ue 


0.1 


0.1 


0.1 




0.8 


0.8 


0.8 


Ymax (electrons) 


1 X 10*^ 


1 X 10^ 


1 X 10^ 


max 


2 


2 


1.5 


Shell width (m) 


1 X 10^ 


1 X 10^ 


1 X 10^ 



4. Discussion 

We saw in the results section that with exactly the same parameters except the jet luminosity 
we get very different lightcurves. This can can be explained by the fact that for the same number 
of shells, increasing the jet kinetic luminosity increases the mass per shell. This in turn has an 
influence on the energy generated at each collision i.e. greater mass shell will have a larger energy 
density. The energy density is the driving force behind the optical depth for a given frequency in the 
shells and thus giving rise to a different spectrum. In short, more massive shells remain optically 
thick to radio frequencies for longer and thus completely changing the overall flux from the jet. 

If we loook at figure ^ we see that having a weaker jet and thus lower energy density for 
each shell leads to most shells becoming optically thin to radio very quickly and thus radio flux 
overall being greater than the IR flux. This can be understood by looking at the position of the self 
absorbed synchrotron spectrum peak relative to the frequency in question 

For the most part T,„ax was kept constant, except in the case of spectrum in figure ^. T,„ax, jet 
kinetic luminosity and shell width have a similar effect on the overall flux. The influence of these 
parameters is almost degenerate in the sense that they influence the energy density of the shells. 
For two colliding shells, greater the difference in their Fs the greater the amount of internal energy 
generated. 

It should be noted that Uth has been set high in the above simulations in order to facilitate the 
longitudinal shell expansion. This expansion should cause the shells to become optically thin to ra- 
dio quicker. We see that the fluxes shown in all the spectra are not in very good agreement with the 
observations for astrophysical systems with similar parameter values. In the case of radio flux es- 
pecially it is highly non-trivial to obtain a flat spectrum with high fluxes. This indeed an interesting 
problem and as noted earlier difficult to solve without making certain assumption either about the 
re-acceleration mechanism [|I|] or the jet geometry [^. In our model, with many collisions taking 
place all along the jet, the shells are continually compressed and powerlaw distributions completely 
replaced thus becoming optically thick to radio. Shells then have to go through adiabatic expansion 
in order to become optically thin to radio. However, the adiabatic energy losses mean that the radio 
flux is significantly lower. We are currently exploring the parameter space that may lead to this 
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somewhat fine balance between adiabatic losses and high radio flux for a flat spectram. 

Although it was not shown in these proceedings, the ejection frequency of the shells also plays 
an important role in the lightcurve and spectra obtained from the jet. The ejection frequency set at 
1 Hz seems to give reasonable results and we not that it is approximately the break frequency in 
X-ray power spectra of black holes in bright hard state (e.g. Cygnus X-1). 

The flaring seen in both radio and IR lightcurves is influenced by various parameters. The 
flaring behaviour is determined by the balance between how quickly the shells lose their energy via 
adiabatic losses and the re-energization/compression due to collisions. Rapid adiabatic losses (i.e. 
larger fraction of energy given to Uth) causes greater flaring in the hghtcurve. 

5. Future directions 

As mentioned in the previous section, our initial aim is to explore the parameter space that 
gives rise to the Radio-IR flat spectrum. This is an attempt to break the degeneracy present in 
some of the parameters. Once we have a reasonable picture of what conditions are required to give 
rise a flat spectrum, we can look at how X-ray variability maybe the driving force behind the shell 
ejection variability. We also aim to add other physical process to the model such as synchrotron 
self Compton radiation and radiative losses for the electrons. 
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